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Abstract 


\J 

Based  on  current  knowledge,  an  assessment  of  the  role  of  ions 
in  the  stratospheric  aerosol  formation  process  is  made.  From  our 
order-of-magnitude  estimate  it  is  shown  that  the  impact  of  the 
average  stratospheric  ion  population  on  aerosol  formation  processed 
in  the  unperturbed  stratosphere  is  not  significant. 


Introduction 


Aerosol  can  be  introduced  into  the  stratosphere  (1)  by  direct 
transport  from  the  troposphere  (including  volcanic  injection), 

(2)  as  a result  of  downward  flux  of  meteoritic  material,  or  (3)  by 
gas  to  particle  conversion  processes  occurring  in  the  stratosphere. 

The  most  widely  accepted  explanation  for  the  existence  of  a sulfate 
particle  layer  in  the  stratosphere  ("Juage  layer")  is  that  most  of 
the  sulfate  particles  are  produced  iui  situ  by  the  oxidation  of  sulfur 
dioxide  into  sulfuric  acid  vapor  which  then  undergoes  either  vapor 
transfer  condensation  to  pre-existing  particles  or  a phase  transition 
to  yield  new  sulfuric  acid  droplets.  It  is  believed  today  that  gaseous 
3ulfur  compounds  reach  the  lower  stratosphere  by  both  the  continual 
diffusional  transport  and  by  sporadic  direct  volcanic  injections; 
mainly  sulfur  dioxide  (Junge,  1971*).  However,  Crutzen  (1976)  recently 
has  suggested  that  in  the  absence  of  violent  volcanic  activity  in- 
sufficient 3ulfur  dioxide  reaches  the  stratosphere  directly  from  the 
troposphere  to  account  for  the  3ulfate  layer.  He  has  proposed  that 
carbonyl  sulfide  (CSO)  from  the  troposphere  is  photolyzed  in  the 
stratosphere  to  produce  largely  atomic  sulfur  which  then  reacts  with 
02  to  form  sulfur  dioxide.  The  principal  initiating  reaction  responsible 
for  the  oxidation  of  sulfur  dioxide  to  eventually  sulfuric  acid  in  the 
stratosphere  is  thought  to  be: 

k 

S02  + OH  + M ■*  HS03  ♦ M 
with  k » (2  to  5)  x 10~31  cm6‘sec-1 

(Davis  197U,  Castleman  197*0 . Other  initiating  reactions  involving 
0(3P),  02(ld),  H02,  NH3,  CH302,  direct  photooxidation,  etc.,  have  been 
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dismissed  since  they  do  not  significantly  contribute  to  sulfuric 
acid  formation. 

The  subsequent  steps  HSO-j  -►  have  not  been  identified  and 

established  yet  and,  until  proven  otherwise,  the  initiating  steps 
SO2  ♦ OH  is  assumed  to  be  rate  determining  for  the  formation  of 
sulfuric  acid. 

The  gas  to  particle  conversion  process  (phase  transition)  for 
sulfuric  acid  can  principally  be  initiated  via  two  main  mechanisms: 

I.  Homogeneous,  heteromolecular  nucleation 

a)  as  binary  system  H2SOI*  + H2O 

b)  as  ternary  or  multicomponent  system  such  as 
H2SO4  + HNO3  + h2o 

c)  as  "ion  nucleation" 

II.  Heterogeneous,  heteromolecular  nucleation 

a)  on  the  surface  of  preexisting  "sulfate  particles" 

(Radius  > O.ly) 

b)  on  the  surface  of  preexisting  Aitken  particles 

o 

(Radius  > 50A) 

c)  on  the  "surface"  of  ions  - ion  grovth. 

Further  growth  into  the  particle  size  range  typically  found  in  the 
"Junge  layer"  (0.1  £ R ^ ly)  occurs  via  (heteromolecular)  condensation 
and  coagulation.  The  relative  importance  of  these  two  conversion 
processes  cannot  be  established  with  certainty  at  present  on  the  basis 
of  experimental  data  gathered  in  the  stratosphere  or  deduced  from 
laboratory  experiments . 
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Various  models  have  been  developed  focusing  on  one  or  several 
critical  physical  and/or  chemical  processes  in  an  attempt  to  explain 
features  of  the  stratospheric  sulfate  layer.  For  example: 

a)  One  dimensional  diffusion  models  combined  with  gas  phase 
chemical  kinetics  but  disregarding  phase  transition  processes 
(Harrison  and  Larson  (1971*)  and  Harker  (1975)) 

b)  Homogeneous  and  heterogeneous  nucleation  and  condensation 
models  combined  with  limited  chemical  kinetics  but  disregarding 
transport  processes  (Kiang,  et  al. , (1973),  Stauffer,  et  al. , 

(1973),  Yue,  et  al.,  (1975),  Hofman,  et  al.,  (1975)). 

Both  types  of  models  in  turn  rely  heavily  upon  stratospheric 
chemical  modelling  efforts  predicting  the  concentration  of  trace  con- 
stituents such  as  OH.  Furthermore,  crucial  assumptions  have  to  be  made 
in  both  models  on  the  physical  and  chemical  state  of  the  atmosphere. 

Figure  1 shows  the  results  from  diffusion  models  and  their 
sensitivity  to  changes  of  boundary  conditions  (tropopause  height,  con- 
centration of  sulfur  compounds  near  the  tropopause,  etc.)  and  vertical 
eddy  diffusion  profiles.  Also  included  in  Fig.  1 is  the  approximate 
range  of  stratospheric  sulfate  measurements  as  reported  by  Lazrus 
and  Gandrud  (Spring,  1971).  Certainly  by  changing  those  uncertain 
atmospheric  parameters,  the  model  can  be  "adjusted"  to  agree  with 
experimental  observations.  The  models  predict  steady  state  maximum 
sulfate  concentrations  of  0.08  to  0.27  ug*m-3  (ambient)  or,  more 
accurately,  the  predictions  are  for  sulfuric  acid  concentrations  of 
5-1  x 107  molecules • cm- 3 (at  18  km),  1.8  x 10s  molecules • cm"3  (at  12  km), 

2.2  x 10®  molecules* cm-3  (at  12  km)  and  1.9  x 10®  molecules* cm* 3 (at  17.5  km). 


These  sulfuric  acid  molecules  are  assumed  to  be  converted  immediately 
into  the  condensed  phase  without  specifying  the  conversion  process. 

Some  basic  results  of  model  calculations  dealing  with  gas  to 
particle  conversion  processes  are  shown  in  Fig.  3.  There  are  several 
crucial  assumptions  in  these  model  calculations  which  have  not  yet  been 
fully  verified  by  laboratory  measurements.  They  *1 1 influence  the 
rate  of  formation  of  new  particles.  (For  example,  vapor  pressure  over 
binary  or  ternary  solutions,  surface  tension  of  embryonic  droplets, 
etc.)  However,  order  of  magnitude  estimates  can  be  made,  and,  in 
conjunction  with  sensitivity  tests,  a confidence  limit  for  the  accuracy 
of  predictions  can  be  made  (Kiang,  et  al. , (1975)).  The  general 
conclusion  which  can  be  drawn  from  those  model  predictions  eliminates 
the  "classical"  binary  phase  transition  process  of  HpSO^  + HjO  as 
responsible  mechanism  for  aerosol  formation.  The  critical  sulfuric 
acid  concentration  required  for  homogeneous,  binary  nucleation  would 
be  of  the  order  of  3.7  x 108  molecules • cm-3.  This  concentration  is 
higher  than  the  maximum  concentration  attainable  in  the  stratosphere 
under  steady  3tate  conditions  as  predicted  by  diffusion  models.  Under- 
lying these  models  is  the  assumption  that  the  vapor  pressure  over 
pure  sulfuric  acid  solution  is  10~6  Torr,  (at  25°C)  one  of  the  lowest  re- 
ported in  the  literature. Conditions  for  nucleation  will  be  even  less 
favorable  for  higher  vapor  pressures.  More  likely,  the  I^SO^  molecules 
will  condense  on  preexisting  particles  as  indicated  in  Fig.  3.  Even 
if  the  available  total  surface  area  per  cubic  centimeter  decreases  from 
U.5  x 10"8  cm2*cnr3  (Castleman,  et  al. , (1975))  to  U.5  x 10“8  cm2-cm-3 
(corresponding  to  0.03  ug.m-3),  homogeneous  binary  nucleation  doesn't 


5 


r 

I 

play  a dominant  role  (see  Hamill> et  al.,  1976). 

Most  likely,  however,  the  gas  to  particle  conversion  will  occur  via 
ternary  nucleation  involving  HpSOl  + HNO3  + H2O  (Kiang,  et  al. , 1975). 

This  process  assumes  that  ^-70%  H2SO4  + VlOJl  HNO3  + ^20%  H2O  solution 
droplets  constitute  an  optimum  mixture  for  which  the  total  vapor 
pressure  is  drastically  reduced.  If  the  estimated  thermodynamic  para- 
meters entering  the  ternary  system  (vapor  pressure  and  surface  tension 
data)  are  correct,  then  this  aerosol  formation  mechanism  would  be 
favored  against  any  other  process  put  forth  to  date.  Current  information 
seems  to  indicate  a sufficiently  high  HNO3  and  H2O  concentration  in  the 
stratosphere.  The  critical  concentration  of  H2SO4  required  for  this 
ternary  nucleation  process  is  less  than  105  molecules *cm“ 3 . This  rather 
low  concentration  for  aerosol  initiation  would  warrant  a constant  source 
for  stratospheric  particles.  This  H2SO4  concentration  is  also  by  at 
least  one  order  of  magnitude  lower  than  required  by  Hofman,  et  al.  (1975) 
for  their  proposed  heterogeneous  condensation  model.  Their  model  involves 
Aitken  nuclei  diffusing  from  the  troposphere  into  the  stratosphere  and 
encountering  a gaseous  sulfuric  acid  layer  (106  I^SOij  molecules • cm- 3 max) 
between  16-22  km  altitude.  Condensation  of  H2SO4  on  those  Aitken 
particles  then  leads  to  particle  sizes  and  concentration  profiles  as 
found  in  the  "Junge  layer". 

Some  of  the  pre-existing  stratospheric  aerosol  originates  from 
extraterrestrial  sources.  Newkirk  and  Eddy  (1964)  have  estimated  the 
relative  concentration  of  meteoritic  debris  in  aerosols  (at  20  km)  and 
concluded  that,  over  a size  range  of  0.1  to  2y  radius,  it  represents  less 
than  10  per  cent.  A completely  "particle  free"  (R  > O.ly)  stratosphere 
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is  therefore  an  unreasonable  assumption. 

It  has  become  obvious  that  most  of  our  understanding  of  stratospheric 
aerosol  formation  leans  heavily  on  model  predictions.  Models  have  not 
yet  been  fully  validated  for  the  stratosphere  and  hence  contain  physical 
and/or  chemical  parameters  that  might  be  quite  different  from  reality. 

The  above  attempt  to  summarize  the  current  opinions  on  stratospheric 
aerosol  formation  processes  constitutes  a necessary  prerequisite  for 
critically  assessing  the  role  of  ions  in  this  important  aspect  of 
stratospheric  chemistry.  The  understanding  of  stratospheric  ion 
chemistry  is  incomplete,  mostly  derived  from  laboratory  experiments  and 
ionospheric  research  work,  but  not  yet  validated  by  stratospheric  in  situ 
measurements.  It  is  within  this  framework  of  uncertainty  that  the  role 
of  ions  in  stratospheric  chemistry  (aerosol  formation)  is  being  dis- 
cussed and  compared — on  a competitive  basis — with  homogeneous  and 
and  heterogeneous  reactions  of  the  electrically  neutral  atmosphere. 

Ions 

The  formation  and  nature  of  ions  of  atmospheric  importance  has  been 
summarized  by  Mohnen  (1971,  1971*)  and  Ferguson  (1975).  The  evolution 
and  final  composition  of  atmospheric  ions  strongly  depends  upon  their 
collisions!  interaction  with  the  neutral  atmosphere  and  thus  becomes  a 
function  of  altitude  (and  temperature ) . The  equilibrium  concentration 
of  ions  is  governed  by  a budget  equation: 

^ » q - an+*n”  - n^+  or  S ( R ) * f(R)dR 

dt  R 

min 
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(1) 
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n(+  or  -).  ^on  concentration 
q:  ion  formation  rate 

a:  recombination  coefficient  for  positive  and  negative  ions 

8(R):  ion  annihilation  coefficient  with  atmospheric  aerosol 

f(R):  aerosol  size  distribution  function 

R:  radius  of  aerosol  particles 

Zikmunda  and  Mohnen  (1972)  have  evaluated  the  integral  0(R)*f(R)dR 

“min 

for  the  stratosphere  and  found  it  to  be  much  smaller  than  the  ion  loss 
due  to  recombination  for  the  unperturbed  stratospheric  particle  load. 

The  ion  formation  rate  q,  derived  from  cosmic  radiation,  is  shown  in 
Fig.  1*. 


Typical  values  between  20  and  50  ion  pairs -cm-3 -sec-1  can  be  expected 
within  the  altitude  range  of  10-20  km,  with  a maximum  around  12  km.  The 
recombination  coefficient  is  of  the  order  of  1 x 106  cm3*sec-1  which 
leads  to  average  ion  equilibrium  concentrations  (n  = 

UU00  < (3Lj-15  < 7000  ions -cm-3 


Actual  measurements  in  the  stratosphere  (summarized  by  Hake  (1973)) 
reveal  ion  profiles  that  follow  the  ionization  profile  and  fall  within 
the  concentration  range  of  1000  - max  10,000  ions* cm-3.  For  our  purpose 
of  ion  assessment,  we  will  assume  two  ion  concentration  values: 
n = 1000  cm-3  and  n = 5000  cm~3  and  thus  cover  minimum  and  average 
conditions.  The  average  lifetime  (e_1  value)  of  ions  in  the  stratosphere 
thus  becomes  1000  sec  and  200  sec,  respectively.  As  mentioned  earlier, 
the  chemical  composition  of  ions  depends  upon  the  gas  composition  and 
the  "reactivity"  of  these  molecules  with  particular  ion  species.  The 
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average  stratospheric  concentration  profiles  for  gases,  particles  (large 
particles  and  condensation  nuclei  or  "Aitken  particles")  and  ions  are 
illustrated  in  Fig.  5.  Ions  can  conceivable  alter  their  chemical 
nature  as  they  collide  with  neutral  gas  molecules.  From  this  emerges 
the  concept  of  "initial  ions  -*■  intermediate  ions  -*■  terminal  ions."  The 
boundary  conditions  of  a limited  average  lifetime  for  ions  determines 
the  minimum  concentration  level  for  which  ion-molecule  interactions 
can  occur.  Assuming  an  ion-molecule  reaction  rate  constant  of 
5 x 10~*°  cm3*sec~3  fixes  this  lowest  trace  gas  concentration  to 
3 x 106  molecules • cm-3  (for  1000  sec  average  ion  lifetime)  as  indicated 
in  Fig.  5.  For  a more  detailed  consideration,  measured  reaction  rate 
constants  for  a specific  ion  molecule  process  should  be  used.  Those 
rate  constants  available  to  date  have  been  compiled  by  Ferguson  (1975). 

Our  current  understanding  of  the  nature  of  positive  ions  at  20  km 
altitude  as  derived  from  measured  reaction  rate  constants  and  measured 
and/or  estimated  stratospheric  molecular  concentrations  (Table  l)  is 
that  the  terminal  ion  is  predominantly  of  the  type  with 

n = U (*0.W),  n = 5 (£35*),  n = 6 (%60.6<)  and  n = 7 (4*).  One  can 
expect,  due  to  the  high  collision  frequency,  that  nitrogen  molecules 
are  also  clustered  to  these  terminal  ions.  One  of  the  unresolved 
problems  in  stratospheric  ion  chemistry  relates  to  the  amount  of  gaseous 
ammonia  present  in  the  stratosphere  since  it  reacts  via  proton  transfer 
with  the  hydronium  ion  to  form  mixed  clusters  of  the  type  (tL^O  )Q 

(this  type  of  ion  constitutes  the  predominant  terminal  ion  core  in  the 
lower  troposphere).  If_ the  ammonia  concentration  in  the  stratosphere  is 
below  3 x 106 , then  no  collisional  interaction  occurs  on  the  average. 
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For  negative  ions,  the  transition  from  initial  ions  to  terminal 
ions  is  considerably  more  complicated  than  for  positive  ions  and  involves 
most  every  molecular  specie  listed  in  Fig.  5-  The  simplified  ion 
evolution  sequence  val'.d  in  the  stratosphere  is  currently  thought  to 


be  (Ferguson  (1975)): 


_ lo2) 

Electron  -+■  0^  * 


[H?0] 


'v-lQ  6 sec 


V10“5  sec  f [Oj 


U x 10~6  sec 
[C05] 


'ulOO  sec 


[hno3] 


[clustered  molecules] 

Ion  Nucleation 

Ion  induced  nucleation  is  a well  established  laboratory  phenomena  in 
atmospheric  sciences  although  it  does  not  play  any  major  role  in  cloud 
physics  in  the  lower  part  of  the  troposphere.  Ion  nucleation  occurs 
principally  at  a supersaturation  ratio  lower  than  is  required  for 
homogeneous  nucleation  because  of  the  reduction  in  the  free  energy 


barrier. 
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In  the  stratosphere,  ion  nucleation  can  only  occur  via  heteromolecular 
nucleation,  i.e.,  involving  two  or  more  gaseous  components  such  as 
H^SO^+H^O  or  H^SO^+KNO^+H^O,  etc. , since  the  supersaturation  ratio 
normally  never  reaches  or  exceeds  3.8,  the  required  value  for  a pure 
water  system. 

The  simplified  equation  describing  the  free  energy  for  the  formation 
of  embryos  in  a binary  (ternary)  system  is: 

AG  = -[A]-kT*ln  JiM  - [B]*KT*ln  J£!  - ([C]-kT-ln  d[r^) 

P[A]  P[B]  ^[C] 

(2) 

+ 4wr2*a  + Sf 

2r (1-1/e) 


[A],  [B],  ([C]):  Molecular  number  concentration;  for  example,  number 

of  water, — sulfuric  acid, — (nitric  acid)  molecules 
PtA],  P[BJ.  (P[C]):  Actual  partial  pressure  of  component  A,  B,  (C) 

P*A]»  P*b]»  (PTc]>:  E(luilibrivun  vapor  pressure  for  the  mixture 

Q:  Ion  charge,  normally  one  elementary  charge 

r:  Cluster  radius 

e:  Dielectric  constant 

A general  solution  of  this  ion-induced  heteromolecular  nucleation 
process  as  a function  of  relative  humidity  is  not  available  to  date. 
However,  Castleman  and  Tang  (1972)  have  extensively  studied  ion  nuclea- 
tion in  a pure  water  system  ( unimolecular ) and  some  of  their  general 
conclusions  can  be  applied  here:  Tne  steady  state  concentrations  of  ion 

''lusters  do  not  increase  continuously  from  the  small  hydrates  up  to  and 
including  the  critical  (stable)  embryo  as  the  supersaturation  ratio 
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increases.  The  existence  of  a continuously  increasing  concentration 
is  not  a prerequisite  to  their  participation  in  the  nucleation  process. 
Conceivably,  the  free  energy  barrier  for  nucleation  can  be  different 
for  different  ion  species  irrespective  of  their  sign  (positive  or 
negative).  The  structure  and  bonding  of  the  small  ion  clusters  are 
the  physical  and  chemical  factors  leading  to  the  macroscopically  observable 
different  values  of  supersaturation  ratios  required  for  different  ions. 

It  becomes  obvious  from  the  work  of  Castleman  and  Tang  that  a single 
term,  Q/r-(l-l/e)  cannot  adequately  describe  binary  or  ternary  ion 
nucleation  systems,  and  is  at  the  most  suitable  for  obtaining  order  of 
magnitude  estimates. 

Keeping  these  constraints  in  mind,  Wiendl  (197*+)  has  estimated  the 
AG  values  (from  equation  2)  for  the  clustering  of  HgSO^-I^O  around  ions. 

The  important  conclusion  that  can  be  derived  from  this  study  is  that 
these  cluster  ions  (as  found  by  Castleman  (1971*)  in  the  pure  water  system) 
do  not  increase  continuously  up  to  and  including  the  stable  embryo. 

One  can  also  deduce  from  Wiendl' s results  that  under  stratospheric 
conditions,  the  number  of  HgSO^  molecules  clustered  to  the  (positive 
or  negative)  ions  must  be  less  than  three:  Wiendl  calculated  the 

^SO^-HgO  ion  cluster  distribution  for  25°C,  relative  humidity  503t , 

Ptotal  = Torr,  = 10“11  Torr  (actual  partial  pressure,  the 

saturation  vapor  pressure  over  pure  HgSOk  solution  wa3  assumed  to  be 
IQ-6  Torr)  and  found  the  most  probable  ion  cluster  to  contain  two 
^SOij  molecules  with  U-12  water  molecules  attached  (peaking  at  n^  2o  = 9)- 


Kiang  and  Stauffer  (1973)  (see  also  Hamill,  et  al.,  1976)  have 
obtained  order  of  magnitude  estimates  for  ion  induced  nucleation  rates 
(based  on  equation  (2))  in  the  binary  HpSOij-^O  system  and  found  the 
minimum  H2SO4  concentration  required  for  ion  nucleation  to  be  by  a factor 
of  ten  lower  than  without  ions,  i.e.  , under  stratospheric  conditions 
of  the  order  of  3 x 10  7 molecules • cm-3  (at  1%  relative  humidity,  de- 
creasing with  increasing  rh  (See  text  to  Fig.  3)).  When  calculating 
the  characteristic  time  for  ion-induced  binary  nucleation,  one  must 
keep  in  mind  that  it  depends  on  the  number  of  ions  available.  Once 
the  H2SOI4  concentration  required  for  ion  nucleation  is  reached,  al  1 
ions  will  simultaneously  become  "stable  embryos"  and  grow  into  larger 
particles.  Therefore  the  total  number  of  stable  embryos  is  limited 

to  the  number  of  ions  present.  The  size  of  these  critical  embryos  is 

o o 

in  the  range  of  1*  to  flA.  Assuming  = 5A  leads  to  a specific 

surface  area  of  A = 3 x 10“ 11  cm2 • cm- 3 for  1000  ions •cm'3  and 
1.5  x 10~10  cm2*cm~3  for  5000  ions’cm-3.  This  yields  (see  text  to 
Fig.  3)  for  1000  ions* cm-3 


-*•  x 


ion 


/— Y*  x A-> 

\2'1,*nH2S0lj/ 


and  for  5000  ions*  cm" 


6 x 106  sec 


t = 1.2  X 106  sec.  (See  Fig.  3) 
ion 

The  influence  of  relative  humidity  dependence  on  r ^ enters  through 
the  size  of  the  critical  (stable)  ion  embryo  but  changes  between 
1 < rh  < 10)t  axe  not  expected  to  cause  more  than  one  order  of  magnitude 
change.  Comparing  these  characteristic  times  for  ion-induced  nucleation 
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with  the  characteristic  time  for  the  ^SO^-HgO  ion-free  system  (Fig.  3) 
forces  us  to  the  conclusion  that  ion-induced  nucleation  does  not  play 
a dominant  role  in  stratospheric  nucleation.  Except  for  very  low 
humidities,  it  constitutes  the  least  likely  mechanism  and  yields  orders 
of  magnitude  higher  characteristic  times  than  the  condensation  process 
on  preexisting  particles,  I_f,  under  perturbed  stratospheric  conditions, 
the  ion  density  locally  and/or  regionally  increases  by  more  than  one 
order  of  magnitude,  then  ions  could  compete  with  the  homogeneous 
binary  nucleation.  But  in  order  to  become  the  dominant  phase  transition 
mechanism,  one  would  have  to  additionally  postulate  a very  low 
atmospheric  total  particle  load  (less  than  10"9  cm2-cnr“3). 

Whether  or  not  the  postulated  ternary  nucleation  process  involving 
H2SO1J  ('vT0!t)  - HNO3  (V10£)  - H2O  {^20%)  exists  is  still  uncertain. 

Unlike  the  case  for  the  f^SO^-^O  system,  there  exist  to  date  no 
experimental  data— either  laboratory  simulation  experiments  or  thermo- 
dynamic data  on  vapor  pressures — on  this  ternary  system.  However,  as 
far  as  ion  involvement  is  concerned,  this  unresolved  problem  has  little 
consequence  on  the  conclusions  made  so  far  on  ion-induced  nucleation: 

l)  The  minimum  I^SOi*  concentration  required  for  the  ternary  nuclea- 
tion is  of  the  order  of  10s* cm" 3 under  stratospheric  conditions.  During 
the  average  lifetime  of  an  ion  (200-1000  sec),  it  yould  not  collide 
with  any  molecules,  and,  therefore,  ion  nucleation  in  the 

stratosphere  would  have  to  be  treated  as  a unimolecular  process  involving 
H2O  only  and  requiring  supersaturation  ratios  (with  respect  to  water) 
in  excess  of  3.8.  In  principle,  one  must  postulate  for  any  ion-induced 
nucleation  process  involving  one  or  more  molecular  species  (of  low  volatility) 


in  addition  to  water  that  the  molecular  species  is  (are)  clustered  to 
the  ion.  This  requires  a minimum  concentration  of  at  least 

4 

3 x 106  molecules • cm-3.  Should  the  stratospheric  concentration  of 

be  below  this  level  (as  it  would  be  if  the  ternary  (10%)  - 

HNO^  (1 0%)  - i^O  (20%)  system  exists,  since  the  equilibrium  vapor 
pressure  over  a mixed  solution  droplet  of  this  composition  is 
<10-6  q>orr),  then  ion-induced  nucleation  would  not  occur  at  all. 

2)  The  boundary  condition  for  ion  nucleation  is  that  there  can 
only  be  as  many  stable  embryos  formed  as  there  are  ions  present.  This 
automatically  establishes  and  limits  the  characteristic  time  T^on 
around  >106  seconds  for  ion  densities  normally  present  in  the  un- 
perturbed stratosphere. 

Ion  Growth 

The  existence  of  ion  growth  is  an  experimentally  well  documented 
fact  in  the  lower  part  of  the  troposphere.  Electrical  mobility 
measurements  of  small  ions  seem  to  correlate  well  with  the  "age"  of 
ions ; a longer  interaction  time  with  the  neutral  gaseous  environment 
yields  lower  mobility  values.  This  phenomena  is,  of  course,  easily 
explainable  with  the  concept  of  ion  evolution  (Fig.  5),  i.e.,  collisional 
interaction  with  various  trace  gases  of  lower  and  lower  concentration 
a3  the  lifetime  of  ion  increases  (Mohnen  (197*0 ).  Ttypical  mobility 
measurements  in  the  lower  troposphere  reveal  a broad,  continuous 
"spectrum"  peaking  around  1.2  cm2-v~**s”3.  The  initial  ion  mobility  is 
usually  in  excess  of  2 cm2 ‘s-1 * v-1 . While  ion  mobility  measurements  do 
not  uniquely  relate  to  the  mass  of  ions,  they  nevertheless  are  a good 
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indicator  for  obtaining  some  information  on  the  average  ion  mass. 

Fig.  6 shows  these  ion  nobility  measurements  for  the  stratosphere, 
ranging  from  0.5  to  1.5  cm2 ■ v-1 • s~ 1 at  20  km.  Whether  or  not  this 
spread  is  real  or  exhibits  artifacts  of  individual  instruments  is  subject 
to  some  debate.  There  are  only  a very  few  mobility  data  available  from 
the  stratosphere.  If  we  assume  that  the  terminal  stratospheric  positive 
ion  is  of  the  type  H^O*- (H20)n , with  n = 6 and  7 (at  20  km),  then  this 
would  yield  a single  mobility  value  of  around  1.75  cm2 • v“ 1 • s“ 1 for 
S.T.P.  (Mohnen  (197*0). 

A departure  towards  lower  values  would  be  indicative  of  molecules 
clustered  to  this  core  ion.  The  most  likely  candidate  is,  of  course, 

N2.  (Although  its  clustering  ability  is  very  low,  the  high  collision 
frequency  would  still  allow  for  some  attachment.)  We  are  again  forced 
to  retreat  in  the  area  of  simple  model  calculations  based  on  laboratory 
measurements  in  an  attempt  to  assess  the  existence  and  importance  of  ion 
growth  in  the  stratosphere.  To  date  there  have  been  two  mechanisms 
observed  by  which  ions  can  grow: 

1)  The  usual  ion  molecule  cluster  development  whereby  the  cluster 
size  and  the  nature  of  clustered  molecules  reflect  and  are  in  a pseudo 
equilibrium  with  the  surrounding  gaseous  environment.  The  forces  by 
which  every  molecule  is  held  together  within  the  cluster  are  hydrogen 
bonding  and/or  electrostatic.  Their  growth  is  limited  by  the  free  energy 
barrier  for  nucleation. 

2)  Chemical  reactions  that  occur  between  clustered  molecules,  or 
chemical  reactions  between  ion  clusters  and  impinging  gas  molecules  whereby 
the  reactant  product  is  incorporated  within  the  ion  cluster.  There  is 
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virtually  no  information  available  on  this  type  of  ion  chemistry  except 
Kadlecek's  laboratory  experiments  where  ion  growth  was  established  in  a 
system  containing  air  and  variable  amounts  of  water  vapor,  ammonia  and 
sulfur  dioxide.  The  ions  that  were  mass  identified  and  mobility  analyzed 
were  in  general 

NH£  • (NH3)n  • (H20)m  • n,m  = function  of  [NH-j],  [H^O ] 

and  NH£  • (NH3)i  • (HgOjj  • [(NH3)2  • S02]k 
Individual  mobility  peaks  were  observed,  ranging  from 
2.1*6  cm2*v“1,s“l  (for  NH£  *(NH3)2  3)  down  to 

1.17  cm2  * if" 1 • s" 1 (for  NH^  * (NH3)3  u)  * l(BH3)2  • S02]U  and  lower. 

o 

At  the  same  time,  some  evidence  of  particle  formation  of  size  £ 25A 
radius  was  obtained,  the  total  number  of  particles  being,  however,  by  at 
least  two  orders  of  magnitude  lower  than  the  ion  concentration.  For 
these  types  of  ion  molecule  reactions , there  appears  to  exist  no  energy 
barrier  to  cluster  growth.  As  these  "crystal-like  structured"  clusters 
develop,  the  addition  of  more  mass  enhances  the  stability  of  the  cluster 
(Castleman  (1973)).  McLaren,  et  al. , (1973)  have  determined  reaction 
enthalpies  for  (NH^g’SOg  of  -1*5  Kcal-mole-1  and  -3**  Kcal*mole_1  in  the 
presence  and  absence  of  water  vapor. 

In  a system  containing  air,  ammonia  and  hydrochloric  acid,  Mohnen 
and  Coffey  (1972)  have  found  particles  to  be  produced  at  ammonia  concen- 
trations as  low  as  100  parts  per  trillion.  In  the  absence  of  ions, 
particle  formation  ceased.  They  postulated  an  overall  ion  reaction 
scheme  NIf£*  (NH3)i*  (ILjO ) j + n*HCl  -*■  particles  as  being  responsible  for 
the  described  phenomena.  Fehsenfeld  (personal  communication,  1975)  has 
found,  however,  that  the  postulated  ion  molecule  reaction  is  endothermic 
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at  least  for  all  i<3  and  J<2.  The  exact  mechanism  of  this  ion-induced 
aerosol  formation  process  is  therefore  not  known  at  this  time.  Growth 
processes  involving  negative  ions  and  leading  to  particle  formation  have 
not  been  reported  to  date,  which  does,  however,  not  exclude  their  existence. 
There  are  conceivably  other  atmospheric  trace  gases,  particularly  reactive 
hydrocarbons , ozone  and  oxides  of  nitrogen,  that  could  undergo  a chemical 
transformation  after  attachment  to  the  ion  cluster  and  cause  it  to 
grow  slowly  into  a larger  and  more  stable  complex.  The  conditions  sire 
favorable  for  such  processes  in  the  lower  part  of  the  troposphere.  (Most 
of  the  trace  gases  mentioned  above  are  present  in  concentrations  ex- 
ceeding 1010  molecules -cm- 3 and  hence  allow  for  at  least  103  collisional 
interactions  during  the  average  lifetime  of  ions.)  But  not  even  there 
does  the  average  ion  grow  beyond  a mass  corresponding  to  *'•1.2  cm2-v-1*s-1. 
Ions  having  a mobility  lower  than  0.5  cm2-v~1sec~1  constitute  less  than 
5%  of  the  total  ion  population. 

In  the  stratosphere  the  concentration  levels  of  trace  gases  such  as 
NH3,  SO2,  reactive  hydrocarbons,  etc.,  are  considerably  lower.  As  mentioned 
above  in  the  section  on  "Ions,"  the  average  positive  and  negative  ion 
will  not  collisionally  interact  with  any  molecular  constituents  having  a 
concentration  less  than  3 x 10s  molecules* cm-3.  Although  the  lfH3  concen- 
trations in  the  stratosphere  have  not  been  measured  yet,  current  estimates 
from  various  investigators  (Crutzen,  personal  communication,  1975)  tend 
to  be  values  less  than  5 x 107  molecules • cm-3  (20  km).  If  this  is  the 
case,  then  ion  growth  due  to  the  (Nl^g'SOg  complex  (the  only  mass 


spectrometrically  identified  ion  growth  process  reported  to  date)  is 
very  unlikely  to  occur  in  the  stratosphere  since  it  requires  first  the 
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formation  of  NHf  (NH_)  core  ions.  This  is  only  achievable  if 
4 on 

> 108  molecules • cm-3.  For  the  same  reason,  one  can  exclude  the 
ion-induced  NH^-HCl  particle  formation  reaction.  Any  chemical  reaction 
between  OH  and  ions  is  also  unlikely  to  occur.  (See  Table  1 and  Fig.  5 
for  the  concentration  range  of  OH. ) 

One  can  conclude  that  ion  growth  leading  to  larger,  more  complex 
and  stable  ions , and  eventually  to  particle  formation , is  not  favored  in 
the  unperturbed  stratosphere. 

Even  if  there  were  yet  undiscovered  chemical  reaction  mechanisms  by 
which  ion3  could  grow  into  particles  that  are  stable  after  the  recombination 
process,  then  their  final  fate  would  be  coagulation  with  preexisting 
particles  rather  than  condensational  growth  into  large  "sulfate  particles": 

T = 1 

coagulation  K 2.n1 

2:  coagulation  coefficient  £ 10-6  cm3* sec”1  for  Rg  £ 5A 

preexisting  particles  I 1 cm-3  for  * O.lu 

Tcoagulation‘  ^ ^ 3ec 

This  characteristic  time  has  to  be  compared  with  the  characteristic  time 
for  condensational  growth. 

tc  = (- — ^ x A'1  = 2.3  x 108  sec 

o 

for  A = 2.5  x 10~13  cm2'cm~3  (5A  Radius  and  100  Particles • cm-3 ) . 

We  assume  here  that  only  a fraction  of  the  total  ion  population  can  grow 
into  stable  particles  (100  cm-3).  Ion  growth  mechanisms  can  therefore 
not  be  responsible  for  the  formation  of  the  "Junge  layer." 
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Conclusion 

This  ion  assessment  with  respect  to  aerosol  formation  in  the 
stratosphere  is  baaed  on  numerous  assumptions,  model  predictions,  de- 
ductions from  processes  occurring  in  other  parts  of  the  atmosphere,  and 
last  but  not  least,  on  rough  order  of  magnitude  estimates.  With  those 
limitations  in  mind,  one  may  conclude  that  the  impact  of  the  average 
stratospheric  ion  population  on  aerosol  formation  processes  in  the  un- 
perturbed stratosphere  is  not  significant. 

In  order  to  have  a better  understanding  of  the  role  of  ions  in  the 
stratospheric  aerosol  formation  processes,  in  situ  measurements  of  con- 
centrations of  ions,  OH,  SO2,  NH3  and  aerosol  size  and  number  distribution 
must  be  made.  Also,  laboratory  investigation  and  theoretical  study 
on  the  thermodynamic  and  kinetic  properties  of  the  formation  of  ion 


cluster  involving  multi-component  system  are  needed. 
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Fig.  1. 

Marker  (1975):  Chemistry  of  Che  stratosphere  has  reached  steady  state. 

Sulfur  dioxide  mixing  ratio  fixed  at  0.2  ppbv  and 
CS023/CS0a  ] “ 3.3  at  the  tropopause.  Three  different 
vertical  Kg  profiles  are  used  (all  within  the  range  as 
indicated  in  Fig.  2). 

Case  1:  tropopause  fixed  at  15  km  (Kg  from  Uofay  and 

McElroy) 

Casa  2:  tropopause  fixed  at  10  km  (Kg  from  Chang) 

Casa  3:  tropopause  fixed  at  10  km  (Kg  from  Bays  and 
Ollvaro) 

Reference:  Alan  B.  Barker,  J.G.R.  80,  No. 24,  3399-3401  (1975) 

Harrison  and  Larson  (1974): 

Sulfur  dioxide  mixing  ratio  fixed  at  1 ppbv  and  [S0^“]  • 0 

at  the  tropopause.  Tropopause  fixed  at  15  km.  Kg  constant 

2-1  fi 

(5000  cm  *s  ) above  tropopause.  Require  3x10  molecules* 

cm-^s-^  upward  flux  of  SO^  through  tropopause. 

Reference:  H.  Harrison  and  T.  Larson,  J.G.R.  79,  No. 21, 

3095-3097  (1974) 

Both  models  imply  a gas  to  particle  conversion  process: 

/m\  2 7 ^ 

S02  + OH  HS03  ■*  H2SO4  (gas)  -V  SO4  (aerosol) 

They  do,  however,  not  specify  the  aerosol  formation  process 

itself. 

Laxrus  and  Gandrud  (1*973): 

Aircraft  measurements  of  sulfate  concentrations  for  Spring 


1971 
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Cont'd  Fig.  1 

Bars  indicate  the  approximate  range  of  S0^~  at  altitude 
15,  17,  and  19  km.  (25°-65°  latitude.) 

Reference:  A.L.  Lazrus  and  B.U.  Gandrud,  "Progress  report 

on  systematic  study  of  stratospheric  aerosol." 
NCAR,  Boulder,  Colo.  (1973) 
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Figure  2 
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Fig.  3.  Character  1st  1 : times  for  binary  and  ternary  nuclaation  and 
for  heterogeneoua  condensation  on  preexisting  particles  as 
function  of  relative  humidity.  (Temperature  -55°C, 
altitude  £20km) 

Curve  H2SO4  - H2O: 

_2 

Equilibrium  water  vapor  pressure  3x10  Torr,  actual  water  vapor 

partial  pressure  3.5  ppmv. 

Formation  rate  for  H2SO4  from  SO2: 

- iLSOzl  - k*  [S02][0H][r]  - dCH2S043 

dt 

[SO2]  ^0.2  ppbv;  [oh]  - 6xl05ca-3;  k “ 3xl0”31  cm® s'1; 

[M]  - 1.85X103,8  cm-3.  Therefore,  production  rate  for  H2SO4: 
eg  ^ 125  molecules* cm-3. 

Composition  of  stable  embryonic  droplet  n4J0Z  H2SO4,  n«20Z  H2O. 
Critical  size  of  stable  embryo:  6-7  X.  Vapor  pressure  of 

H^SO^  over  pure  solution  10~®  Torr.  H2SO4  gas  phase  concentration 
required  for  the  binary  nucleation  process  to  occur: 
v3.7xlO  molecules*  cm  (at  1Z  relative  humidity,  decreasing  with 
increasing  relative  humidity). 

For  further  details,  see  "The  Interaction  of  Atmospheric  Sulfur 
Compounds  with  Cloud  and  Precipitation  Elements,"  Interim  Report  to 
NSF.  ASRC  Publ.  No.  365,  May  1975.  Authors:  C.K.  Turn,  C.S.  Klang, 

V.A.  Mohnen  and  E.  Danlelsen. 

H2SO4-HNO3-H2O: 

All  parameters  the  same  as  in  the  binary  system.  Critical  composi- 
tion of  stable  embryonic  droplet:  75Z  H2SO4,  10Z  HHO3,  20Z  H20. 


Conc’d  Fig.  3. 
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Under  Chose  conditions,  Che  partial  vapor  pressure  for  chla 
ternary  system  is  estimated  to  be  10-^,  10~^  and  10"7  Torr 
respectively.  The  minimum  concentration  of  f^SO^  required  for  this 
ternary  homogeneous  nucleatlon  process  is  of  the  order  of  2x10^ 
molecules *cm~^  (for  1Z  relative  humidity).  For  further  details  see 
"The  Interaction  of  Atmospheric  Sulfur  Compounds  with  Cloud  and 
Precipitation  Elements"  cited  above. 

PARTICLES; 

It  is  assumed  chat  the  preexisting  "Junge  layer"  contains  0.3  pg/M~^ 
of  particulate  matter  with  an  average  size  of  R % O.lu,  a particle  density 
of  p - 2 g*ca~^  and  a total  number  concentration  of  around  40.  This 
leads  to  a total  surface  area  per  unit  volume  of  A - 4.5x10  cm  *cm~J. 
Deposition  of  H^SO^  molecules  on  these  preexisting  particles 
("condensation")  is  proportional  to  the  surface  area  (this  condition 
is  fulfilled  for  all  particles  with  R<0.3p).  Sticking  coefficient  is 
assumed  to  be  one.  The  characteristic  time  is  then  defined  by 

T particle  % ‘ ^ 

K:  Boltzman  constant,  T absolute  temperature 

m:  mass  of  impinging  gas  molecule,  here  V°4 

^ ” -? ‘2L  x 10~^  cm^*cm-^ 

R‘P 

_3 

y:  particulate  matter  concentration  in  pg*m 

_3 

p:  density  in  g’em  and  R:  particle  radius  in  microns 

(Note  that  Tpart^cle  ~ A-1,  hence,  a 10  fold  decrease  in  specific 
surface  area  will  result  in  a tenfold  Increase  in  characteristic  time.) 
(For  further  information,  see  above  cited  report.) 
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Cont'd  Fig.  3. 

Definition  of  characteristic  time: 

Minimum  concentration  required  for  the  homogeneous 

„ nucleatlon  process  — 

nucleatlon  Production  rate  of  H2SO4 


ALTITUDE  KM 


Figure  4 


o £ 

Q.  Ui 

v>  ^ 

S * 

5 ^ 

*»  * 


\ 


1958  1965 


A 1958  75 * LAT. 

q 1958  55° 

C 1565  55* 

0 1965  75° 


ION  PAIR  PRODUCTION  RATE  CM  * SEC  ' 


Source : H.V.  Neher,  J.G.R.  76,  1637-1651  (1971) 

R.D.  Hake,  E.T.  Pierce  end  W.  Vlezzee,  "Stratospheric  Electricity. 
Final  Report,  SRI  Project  1724,  Menlo  Park,  Ca.  (1973) 

M. A.  Ruderman  and  J.V.  Chamberlain,  Doc. Ho.  JSS-73-18-3,  Institute 
for  Defense  Analysis  (1973) 
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Tab la  1 

AVERAGE  CONDITIONS  IN  STRATOSPHERE  (20  ka) 


Trace  Caa 

Trace  Caa 

Concentration  la 
Molecules *cn_^ 

Characteristic  Ion^ 
Interaction  Tine 
In  Seconds 

No. of  Colllelonal  Interactions 
Durlnx  Average  Lifetime  of  Ion(9 
For  For 

1000  Iona/cir*  S000  Ions/cn 

o2<l> 

3.9xl017 

5xl0"9 

2x1 011 

4,1010 

co2(1) 

5.9X101* 

2.9x10“® 

3.4x10® 

6.8xl07 

10l3 

2xl0~* 

5x10® 

105 

o3  (i) 

2.5x10s2 

8xl0~* 

1.3x10® 

2.6x10s 

N02(l)  ; 

1.5xl09 

0.73 

1.3x10s 

2.6xl02 

NO  (^ 

5x10* 

4 

250 

50 

HNO^1* 

5xl09 

0.4 

2.5x10s 

500 

n2o5 

4x1 06 

500 

2 

- 

no3 

105 

2x10* 

- 

- 

2xl09 

1 

1000 

200 

CIO 

2xl07 

100 

10 

2 

OH 

6xl05 

3.3x10s 

(0.3) 

- 

o2h 

2xl07 

100 

10 

2 

0(3P) 

3x10s 

6.7xl03 

(0.15) 

- 

ch3o2 

lxlO7 

200 

5 

1 

Unknown 

• 

• 

so2(2) 

<3.7x10® 

<3 

<200 

<40 

h2soa(3) 

5 8 

>10  but  <10 

<20 

<50 

<10 

nh3  W 

«3.109 

«1 

«1000 

«200 

Particles^ 
(40  cw“3 
R K 0.1m) 

2.5x10s 

— 

~ 

j 

3U 
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(5)  The  lifetime  of  lone  la  defined  by 

do*  ■ _ u • n~ 
de 

a:  recombination  coefficient,  here  aeauned  to  be  of  the  order  of 

10”*  cm3 • aec“~ 


n+,  n"  concentration  of  positive  and  negative  Iona  per  cm3  respectively 
t+  (1000  ions/ca3)  %1000  sec 
t+  (5000  lons/cm3)  ^200  sec 

(6)  Characteristic  Interaction  time  Is  defined  by 

A typical  Ion  molecule  rate  constant  of  5xl0“10  cm3 -sec  la  uaed.  It  la 
assumed  that  temporary  or  permanent  attachment  of  the  trace  gas  molecule 
to  the  Ion  occurs,  or  In  some  instances  that  an  ion  molecule  reaction 
occurs . (See  text) 

(7)  Ion  annihilation  time  defined  as  (Zitaaunda  and  Mohnen  (1972)): 

x - 

B-H 

g:  attachment  coefficient  of  Iona  to  particles 

0 - 10"5  cm3  sec"1  for  R ^ O.ly 

-3 

N:  total  number  of  particles.  N - 40  cm 


'trace  gas  concentration]  x rate  constant 


AVERAGE  INTERACTION  TIME  FOR  IONS  ( ION  EVOLUT I ON  J 


LOG, 0 NUMBER  CONCENTRATION  CM 
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Figure  6 

DEDUCED  PROFILES  OF  POSITIVE  ION  MOBILITY 


0 1 2 
REDUCED  MOBILITY  - cm2  V"1  «H 


Source:  Hake  et  al.,  1973 


